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The scaling properties of a phase-ordering system 
with a conserved order parameter are studied. The the- 
ory developed leads to scaling functions satisfying certain 
general properties including the Tomita sum rule. The 
theory also gives good agreement with numerical results 
for the order parameter scaling function in three dimen- 
sions. The values of the associated nonequilibrium decay 
exponents are given by the known lower bounds. 
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I. INTRODUCTION 

In the area of phase-ordering kinetics JlJ there 
are a wide variety of systems which satisfy a form 
of dynamical scaling. In the case of systems with 
a nonconserved order parameter (NCOP) we have a 
simple approximate theoretical model due to Ohta, 
Jasnow and Kawasaki § (OJK) which captures the 
main scaling properties of the associated physical 
systems. In the case of a conserved order parameter 
(COP) the situation is much less satisfactory. The 
difficulty in the COP case is that there are com- 
peting length scales which lead to the necessity of 
treating crossover. This crossover connects up the 
short-scaled-distance nonanalytic domain wall be- 
havior associated with Porod's law |3|-^| and the 
large distance constraints of the conservation law. 
A theory is offered in this paper which is consistent 
with all of the prominent scaling features in the case 
of a conserved scalar order parameter. 

An auxiliary field method is used in essen- 
tially all the available explicit calculations of scal- 
ing functions in phase-ordering kinetics. Thus in 
the OJK approach IjjH-S], and our previous work 
]To| p^| , a local mapping from the original order 
parameter ip onto an approximately gaussian vari- 
able m was developed. While these methods can be 
shown to work well for the NCOP case, there are 



severe limitations in the COP case. In particular 
Yeung, Oono and Shinozaki jl3) found that such a 
local mapping in the COP case led to mathemat- 
ically unacceptable results within the theory. All 
of these theories were developed with the idea that 
the mapping, tp — > ip(m), leads to an equation of 
motion for m which could be argued to be approxi- 
mately consistent with a gaussian distribution for m. 
In Refs.( 0) and ( (ll|) a different approach was 
taken in treating the NCOP case. It was shown that 
the equation of motion for the auxiliary field m can 
be constructed with appeal only to the form of the 
growth law, defined below, and certain general sym- 
metry constraints. Using a novel expansion method 
it was shown how one could obtain the OJK result 
in zeroth order of a systematic expansion. Higher 
order corrections for the associated nonequilibrium 
exponents were also obtained. 

In this paper this general idea is applied to 
the simplest COP system given by the Cahn-Hilliard 
Jl6|] (CH) model. One is led to introduce a nonlocal 
mapping between the order parameter and the auxil- 
iary field and only certain general constraints on the 
theory, like the conservation law and the generalized 
form of the Tomita sum rule [jiTj, are used to de- 
termine the parameters characterizing the nonlocal 
mapping and the correlation function for the under- 
lying auxiliary field. Thus one has a nonlinear selec- 
tion problem where one simultaneously constructs 
the parameters of the mapping and the scaling func- 
tion. The selected scaling function is found to be 
in good agreement with the best available numerical 
results in three dimensions. The nonequilibrium ex- 
ponents are also determined in this theory and take 
on values, as for the OJK theory, corresponding to 
the known (l8) lower bounds. 

II. PHENOMENOLOGY 

The equation of motion in the Cahn-Hilliard 
model governing the conserved scalar order param- 



eter, ip, is given in dimcnsionless form by: 

1 p = v 2 (y'(V')-v 2 v) . 



(i) 



V(ip) is a degenerate double-well potential. Typi- 
cally this system is driven by a set of uncorrelated 
random initial conditions. We look here at critical 
quenches where < tp >= 0. In this case the system is 
unstable and responds by locally growing competing 
degenerate patches of the stable low-temperature- 
ordered phases. These patches correspond to do- 
mains separated by sharp walls of width £. As time 
evolves these domains coarsen and the order grows 
to progressively longer length scales. The growth 
law, L(t), increases without bound with time t after 
the quench. At long enough times L(t) dominates, 
L(t) 3> £, and the order parameter correlation func- 
tion satisfies the scaling equation 19 21 



C(r,t) = (^(r,t)V(0,t)) = ^F{x) 



(2) 



where x = r/L(t), and ipo is the magnitude of %j) in 
the ordered state. The structure factor, the Fourier 
transform of C(r,t), satisfies 



C(q,t) = L*i/>ZF(Q) 



(3) 



where Q = qL is a scaled wavenumber and d the 
number of spatial dimensions. There are a number 
of general properties a correct theory of the phase 
ordering CH model must satisfy [p2||l2| |: 

1. The growth law is given by the Lifshitz- 
Slyozov- Wagner |^3| form: L w t 1 ^ 3 . This result can 
be obtained a number of different ways ]24| that 
are all consistent. 

2. The scaled-structure factor F(Q) has a 
Q 4 behavior f27| for small scaled wavenumbers Q 
directly reflecting the conservation law. 

3. The scaled-structure factor satisfies 
Porod's law 0-0,|2£| for large scaled wavenumbers, 
F{Q) » Q-^K 

4. The scaled-structure factor also satisfies 
the Tomita sum rules. The scaling function 
F(x) has no even terms, except the first, in its ex- 
pansion in powers of x: 



F(x) = l + F lX + F 3 x 3 



(4) 



While there are theories which satisfy some of 
these requirements, there has been none so far which 
satisfies all four. Our goal here is to find the simplest 
theory for the COP case which does satisfy all of 
these properties and leads to an explicit form for 



the scaled structure factor which can be compared 
to the best numerical determinations of F(Q). The 
challenge is to match the small Q 4 behavior for F{Q) 
with the large Q behavior (in three dimensions) 
while preserving the Tomita sum rule. Our approach 
will be similar to that developed in Refs.( H)) and 
( Jl5|) but with some significant differences. 



III. AUXILIARY FIELD MAPPING FOR COP 

As in previous work, we assume in the scaling 
regime that the order parameter can be decomposed 
into an ordering component, which contributes to 
the order parameter structure factor at 0(1), and a 
fluctuating piece which is of higher order in powers 
of 1/L(t): 



ijj = ^ + Q/L 



(5) 



where O is 0(1). Here we assume that the ordering 
component, ip, can be written in the form: 



ij) = a(m + u) 



(6) 



where, as usual, a(m) is the solution of the classical 
intcrfacial equation 



d 2 a 
dm 2 



= V'(a) 



(7) 



with the boundary conditions ]xm m — > ± oc a(m) = 
±-0o- It turns out in this case that it is neces- 
sary to introduce two independent fields m and u. 
We will organize the theory such that m is treated 
as the auxiliary field which, as in the NCOP case, 
governs the short scaled distance physics associated 
with Porod's law and the Tomita sum rule. We will 
assume, as a first approximation, that m is a gaus- 
sian variable driven by an equation of motion of a 
general form compatible with scaling and with coef- 
ficients determined by selection processes described 
in detail below. The quantity O is, in principle, to 
be determined as a function of the more fundamen- 
tal variables m and u. In practice we will not need 
to construct the correlations for the field u explicitly 
or assume that it is a gaussian field. 

If we insert the ansatz given by Eq.(||) into 
the equation of motion, Eq. ([!]) , we obtain 



dt]j 
~dt 



d_ 

at 



(y'$ + o/L) - v 2 (V> + o/l} 



(8) 
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We can use scaling arguments to estimate the contri- 
butions of various terms. We explicitly assume that 
L(t) = Lot 1 / 3 in the scaling regime. The first term 

on the left-hand side of Eq.(|), ff , is of 0(L- 3 ) in 
the scaling regime. Then we can estimate 



0_ 

dt 



(9) 



and this term can be dropped when compared to 
the leading order in the equation of motion. Next 
we can expand 

v'{4> + e/L) = v'(^) + v"(^)y + o(L- 2 ) . (io) 

In comparing these two terms we have, using Eq.(|^) 
and m ~ L, that 



d 2 a(m + u) mtT-i\ 



dm? 



0{L- 



(11) 



and the term proportional to V"(ip) = V"(tpo) + 
dominates at leading order for long times. 
Finally we have that the last term on the right hand 
side of Eq.(|l|), 



7 4 (vi + e/x) «o(l- 4 ) 



(12) 



and can be dropped in Eq.(g). The equation of mo- 
tion then reduces to the key result 



dip 
~dt 



^v 2 e 



(13) 



where k = V"(ip ) > 0. 

The quantity Q is arbitrary except for the 
very important constraints that it be of 0(1) in the 
scaling regime and it must be consistent with the 
system ordering. It is at this stage that we real- 
ize that there is an apparent flexibility in the con- 
struction of the scaling solution. Since there is a 
belief that the scaling functions are universal, there 
must be mechanisms, like the nonlinear eigenvalue 
problem encountered in Ref . ( |l(J , which selects the 
scaling structures which do not directly depend on 
the physics at the smaller length scales. The key 
assumption in going forward is that the auxiliary 
field method can be used to describe the inner scal- 
ing regime and there is a crossover to an outer scal- 
ing regime dominated by the conservation law. The 
building blocks we can use to construct 0, which are 
compatible with this crossover and are of 0(1), are 
ip and a(m). 



The simplest assumption p9] for the Fourier 
transform Q q (t) , robust enough to give a satisfactory 
solution to the problem, is of the form: 



Q q (t) = M q (t)^ q (t) - N q (t)a q (t) 



(14) 



where M q (t) and N q (t) are functions to be deter- 
mined The equation of motion, Eq.([l3|), is given then 
by: 



dt 



= -M q (t)ij q (t) + N q {t)a q (t) (15) 



where 



MS) = JTqMS) 



NS) = ^NS) 



(16) 



(17) 



We can easily obtain a partial solution to Eq.(15) 
Let us define the auxiliary quantity 



.J2*rAf,(r)) 



(18) 



U q (ti,t 2 ) = e v 
and write 

$S) = u S> t o)xS) ■ (1 9 ) 

Taking the time derivative of this expression we find 



d ^ {t) --M q (t)Mt) + ust ) dXq{t) 



dt 



dt 



(20) 



Comparing Eqs.(|l5[) and fl20|), we obtain the equa- 
tion for Xq(t), 



JT (t t \ dx ^ 

USM-gf- 



ns)*S) 



which can be rewritten as 
dxS) 



dt 



= U q {t a ,t)N q {t)aS) 



(21) 



(22) 



This equation has the solution 



Xq(t) = M*o)+ di U q (t ,t)N q (t)a q (t) (23) 

J to 



ijS) = US,t )i> q (to)+ dt U q (t,t)N q (t)a q (t) . (24) 

Jto 
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This nonlocal relationship between ip q (t) and c q (t) 
should be contrasted with the local mapping used 
in previous theories. Because of the nonlocality the 
criticism of local mappings in the COP case due to 
Yeung, Oono, and Shinozaki [|l3|,[50| is irrelevant for 
the discussion here. 

Notice that we need to determined the func- 
tions M q (t), N q (t), and the variance of the field m. 
Averages over m are discussed in more detail be- 
low. Focusing on N and M, for our purposes here, 
we only need these quantities in the scaling regime. 
Inspection of Eq. (|lj) shows that a general form com- 
patible with scaling is given by 



M q (t) 



dG(Q 2 

at 



where Q = qL{t). Similarly, 

dG (Q 2 ) 



N q (t) = 
We can then write 

8G(Q 2 ) 8Q 2 



dt 



(25) 



(26) 



M q (t) 



dQ 2 dt 



dG{Q 2 ) 8 
dQ 2 dt 



q 2 L 2 t 2 / 3 



where 



Similarly, 



_ dG(Q 2 ) 2 Q 2 = 2 H{Q 2 ) 
dQ 2 3 t ~ 3 t 

H(Q ) - Q —qq^- ■ 
2H (Q 2 ) 



N q (t) 



t 



(27) 



(28) 



(29) 



For our purposes it will be sufficient to assume that 
H and Hq have power series forms: 



H(Q) = J2^Q T 



(30) 



(31) 



For reasons discussed below, we will work explicitly 
with a model where only j 2 , 710, 7 2 , and 7° are 
nonzero. 



IV. STRUCTURE FACTOR 

The quantity of central interest is the order 
parameter structure factor: 

C(q,h,t 2 ) =< tpqfaWSx) > 



(32) 



where in the second line we recognize that in the 
scaling regime only the ordering component of the 
order parameter contributes to the structure factor. 
Inserting Eq.(p3) for ip q (ti) we obtain: 



C(q,tx,t 2 )=< / dti U q (ti,ti)N q (ti)<T 9 (ti) 



x dt 2 U q (t 2 ,i 2 )N q (t 2 )<J q (i 2 ) > 

J t 



dh U q (ti,ti)N q {h) 



dt 2 U q (t 2 , t 2 )N q (t 2 )C (T (q,t 1 ,t 2 ) (33) 



to 



where 



C a (q,ti,t 2 ) =< <J q (h)a q (t 2 ) > 



(34) 



Notice that we have assumed that the correlations 
with the initial state have decayed to zero for long 
times compared to the 0(1) terms contributing to 
the scaling function. We also note that C(q,ti,t 2 ) 
depends only on the magnitude of q. 

The next step is to realize that C a (q,ti,t 2 ), 
which we calculate explicitly below, satisfies a scal- 
ing relation: 



C a (q, ti,t 2 ) = L d {ti,t 2 )iplFa{qL{t x ,h), t u t 2 ) 
Inserting Eq.([35|) into Eq.(f33|) gives 



(35) 



C(q,h,t 2 ) = / dh U q (ti,h)N q (tl) 
'to 



dt 2 [/ g (t 2 ,t 2 )iV g (t 2 )£ d (i 1 ,t 2 )^i ?, CT (gi(ti,«2),«i,t2) • (36) 



Let us now define 
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L 2 T (h,t 2 ) = ^{L 2 (h) + L 2 (t 2 )) 

and choose a time T such that 

L T (h,t 2 ) = L(T) . 
Remembering L(t) = Lot 1 / 3 , we have 

3/2 



T : 



(37) 



(38) 



(39) 



Clearly for ti = t 2 = t, T = t. Now make the change 
of variables t\ — Ts\ and i 2 — Ts 2 in Eq.(j36"l). This 
requires treatment of the quantity 

L T (ii,i 2 ) = L T {Tsi,Ts 2 ) 



x U q {t 2 ,Ts 2 )TN q (T Sl )TN q {Ts 2 ) 



x e d ( Sl ,s 2 )F a (Q£( Sl ,s 2 ),T Sl ,T S2 ) 



(43) 



We can simplify things a bit in the quantities U q , 
defined by Eq.(jl^), where we have in the argument 
of the exponentials 



ft ft/T 

/ dr M q (r) = / dz x M q (T Zl ) 

JTsx J s t 



(44) 



Inserting our general form for M q , given by Eq.(27), 
we have 



Si 



L(T)£(s u s 2 ) = L T (t 1 ,t2)£(s 1 ,s 2 ) 



(40) 



where 



e(si,s 2 ) = \/l(s 3 1 /2 + sl /2 ) . (41) 



Eq.(|36[) then takes the form: 
rti/T 

C(q,t!,t 2 ) = / Td Sl U q {t 1 ,Ts 1 )N q {Ts 1 ) 

Jta/T 
rt 2 /T 

x / Tds 2 U q {t 2 ,i 2 )N q (Ts 2 ) 

Jta/T 

x L d (T Sl ,Ts 2 )F a (qL T {t u t 2 )£( Sl , s 2 ) 7 T Sl ,Ts 2 ) 



rti/T rt 2 /T 

ipl I dsi U q (ti,Tsi) / ds 2 U q (t 2 ,Ts 2 ) 

Jta/T Jt /T 



x TN q (T Sl )TN q (T S2 )L d (T)£ d { Sl ,s 2 ) 



x F a {Ql{ Sl ,s 2 ),Tsx,Ts 2 ) 



= L d (T)^F(Q,t 1; t 2 ) 



(42) 



where Q = qLT(ti,t 2 ) = qL(T) and the physical 
scaling function is given by 



fti/T pt 2 /T 

F{Q,t 1 ,t 2 ) = dsxU q {t x ,Tsi) ds 2 

Jta/T Jta/T 



t 2 /T 
to/T 



2 r t/T dzi 

2 / ^H(Q£( Zl )) 

3 J Sl z\ 



where £{z\) = z^ 3 , and similarly 



TN q (T Sl ) = T- 



2H (Q£( Sl )) _2H (Q£( Sl )) 



Then 



f_ 2 rt/T 

U q {t u T Sl ) = e\ 3j °i n 



H(Q£( Zl ))) 



(45) 



(46) 



(47) 



and the scaling function is given by 

fti/T pt 2 /T 

F{Q,t 1 ,t 2 ) = d Sl U q {h,T Sl ) ds 2 U q {t 2 ,Ts 2 ) 

Jta/T Jta/T 



2 H (Q£( Sl )) 2 H Q (Q£(s 2 )) 



S2 



x£ d ( Sl ,s 2 )F a (Q£( Sl , S2 ),T Sl ,T S2 ) 



(48) 



We can then make one final change of integration 
variables to y-i = s 2 J 3 (i — 1,2).. We have then 



£(si,s 2 ) 



Vi + V2 



' 2/3 . 2/3 

s/ +s 2 < 



1/2 



1/2 



(49) 



the equation for the scaling function becomes, 
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F(Q,h,t 2 )= ( {tl ' T) ^ U q {ti,T y \' 2 ) 

J(t /T) 2 /3 Vl 



< / ^ C/g (t 2 ,T 2/2 3/2 ) 

'(t /T) 2 /3 J/2 



x H (Q^)H (QVm) 
x ^( 2 / 1 ,y 2 )^(g^( 2 / 1 ,y 2 ),Ty? /2 ,r 2/2 3/2 ), (50) 

2/3 

and, with 7/1 — z x in the integral, 



3 \ J V1 



- H(QVW) ) . (51) 



U q (h,Tyl /2 ) = e 



For our simple polynomial model given by Eq.(|3C 
we can carry out the y\ integration in U q explicitly. 



for t > to. Since the growth law is of the form 
L = Lot 1 ' 3 , we see that for a scaling result we 
must choose uj q {t) w O (L~ 3 ). Since we can esti- 
mate q « O (L^ 1 ) in the scaling regime, we write 
quite generally that 



L 2 



= jtt (a Q + caQ + a 2 Q 2 + a 2 Q 2 + a 3 Q 3 ) (57) 

where Q = qL. We can use general arguments, sim- 
ilar to those due to Ohta and Nozaki pl| , to justify 
including the odd terms in Q. Eq.(|57|) can be rewrit- 
ten in the convenient form 



, , dCl(Q(t)) d , 



(58) 



where 



V. AUXILIARY FIELD SCALING FUNCTION 

We must now work out the scaling properties 
of the field a(m(r, t)) in the case where m is a gaus- 
sian variable. It is well known that in the scaling 
regime 

C <J (r 1 ,t 1 ,Y 2 t 2 ) =< a{m(r 1 ,ti))a{m{Y 2 ,t 2 )) > 



ip 2 -sin 1 /o(ri,ti,r 2 ,t 2 ) 

7T 



where 



/o(ri,ti,r 2 ,t 2 ) 



Co(ri,fi,r2,f 2 ) 
A/5o(*i)So(i2) ' 



C (ri,ti,r 2 ,£ 2 ) =< m(ri,ti)m(r2,t 2 ) > 
and 

Sb(*i) = C (ri,ti,ri,ti) . 



(52) 



(53) 



(54) 



(55) 



Thus we need to focus on the determination of the 
auxiliary field correlation function Cq. 

Let us assume that the general equation of 
motion satisfied at gaussian level by the auxiliary 
field m has the local form in Fourier space: 



dm q (t) 
dt 



-ui q (t)m q (t) 



(56) 



fi(Q) =aQ + bQ 2 + cQ 3 



(59) 



= -a Q /L% a = ai/Ll, b = a 2 /(2i§), and c = 
a 3 /(3Lg). 

The equation of motion for the auxiliary field, 
Eq.(pq), has the general solution 



m q (t) = e 



Jt ° m q {t ) 



(60) 



Inserting Eq. ( p8[ ) into the integral in the exponential 
gives 

m (t) = e -(nCQ(i))-nCQ(io)) e ^o J»(*/*o) m 



-(O(Q(t))-O(Q(t )) 



m q (t ) 



(61) 



If we then form the auxiliary field correlation func- 
tion, as defined by the Fourier transform of Eq.(54) 
we obtain 



C (q,ti,t 2 ) = ( 



-(n(Q(ti))+fi(o(t 2 ))-2n(Q(t )) 



xC (g,* ,*o) ■ (62) 
If we go to the equal-time limit this becomes 

L) e -(2nWW)-anW(to)) Co(9>to>to) 
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The important autocorrelation function, defined by 
Eq.(B5h, is given then by 



S (t) 



d d q 



(2tt) 



^Cofa,*,*) 



L d {t) J (2n) d \t 



2u 

— I e' 



-(2O(Q)-2O(QL(t )/L(t))) 



X C (Q/L(t),t ,t ) 

For large L(t) we can replace 

C (Q/L(t),t ,t ) -» g 



(64) 



(65) 



where go is characteristic of the initial correlation 
function, and 



2fi(Q) - 2Q,{QL(to) / L{t)) 2Q(Q) 



(66) 



to obtain in Eq.(p4J) 



So(*) 



,9o 



2u>o 



d d Q 

(27T) rf< 



-(20(Q) 



Here we have introduced the constant 
J 4 = f . 



(67) 



(68) 



By assumption the auxiliary field scales with the 
growth law, vn k, L and So(t) w L 2 (t). Using this 
result back in Eq.(p7n gives 



S (*) « J L 2wo3 - d m L 2 . (69) 
This fixes the constant loq to have the value: 
2 + d 



w 



6 



(70) 



Using this result in the scaling regime, the auxiliary 
field correlation function, given by Eq.(|6^) can be 
written in the form 



Co(q,h,t 2 ) 



t 2 

i n 



Co(g,t 1 ,t 2 ) 
VSoJhWoJh) 



(72) 



Using previous results this can easily be put into the 
form 

fo{qMM) = {L{ti)L{t 2 )f 2 . (73) 

Jd 

With Q,{Q) given by Eq.(p9|), we have in the argu- 
ment of the exponential: 

n(Q(*i)) + n(Q(* a )) = o?(L(ti) + L(i a )) 



+ 6<^(Z%) + L 2 (i 2 )) + ^(* 2 )) . (74) 

If we introduce L T (ti,t 2 ), as in Eq. (pS7|) , and Q 
(j 1. 1 . we obtain 

n(Q(ti)) + n(Q(* 2 )) = 2aQh(t 1 ,t 2 ) 



+ 2bQ 2 e 2 (t ll t 2 ) + 2cQ 3 £ 3 {t 1 ,t 2 ) 



(75) 



where 



1 (L n {t 1 )+L n {t 2 ) 



2 \ LJ^(ti, t 2 ) 



1 /3 + 



(76) 



For equal times ^„(t, i) = 1, while for £i ^> i 2 

4(*i,*2)=2"/ 2 - 1 (77) 

and £ 2 (ti,t 2 ) = 1. We can then choose Lo such that 
c = 1/4, 26 = /i and 2a = —A. We then have the 
final result 



/o(«,ti,* 2 ) = (£(ti)i(*2)) 



d/2 e 



-JV(Q,ti,t 2 ) 



•/,J 



(78) 



where 



N(Q,ti,t 2 



J 4»(ti,*2) + pQ 2 



.(*i,*2) • (79) 



-(Q(Q(ti))+ncQ(i2))) 



.90 ■ (71) 



The inverse Fourier transform of fo{q,ti,t 2 ) is the 
quantity which is related to the a- correlation func- 
tion by Eq.(H): 



Eventually we need the inverse Fourier transform of 
the normalized correlation function 



dd i ~-m-? t 

(27T) 



fo(r,h,h)= j 7 ^-e-^ r fo{q,h,t 2 ) 
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d d q 
(2n) d 



\L( tl )L(t 2 )) 



d/2 e" 



-N(Q,t!,t 2 ) 



J,, 



1 ^ - r ^h(qvw) 



1)2 



x H (Q^y{)H {Q^yi) 



d/2 



L(h)L(t 2 ) 
L 2 T (t,t 2 ) J J (2ir) d ' 



where x = rj Lt and now 



d d Q „- t Q. S e- N(QM) 



■hi 



(80) 



Jd = 



d d C 



d d Q _ 
{2ir) de 



N{Q,t,t) 



(81) 



The Fourier transform of the tr-correlation function 
is given by: 

C a (q,t u t 2 ) = J d d re"** -sin- 1 (/o(r,*i,t 2 )) 



d* x e +i & s -sin- 1 (f (x,t u t 2 )) (82) 

7T 



and the scaling function appearing in our physical 
scaling function, Eq.(|5^), is given by 



£ d [3H,V2)F a (Q£( V1 ,y 2 ),ty 3 1 /2 ,ty 3/ ^ 



(85) 



At first sight, since H(Q^/yi) goes at least as 
fast as Q 2 for large Q, this integral looks like 
it gives exponential behavior in Q for large Q. 
Closer inspection shows the asymptotic depen- 
dence on Q is much slower in those regions of 
the y± and y 2 integrals near 1. By expanding 
£ d (yi,y2)F a (Q£(y 1 ,y 2 ),tyl /2 ,tyl /2 ) about y x = 1 
and y 2 = 1, it is not difficult to show, assuming that 
Hq and H increase algebraically with Q for large Q, 
that 



F(Q) = + O (H(Q)r 1 ) 



(86) 



F tT (Q,t 1 ,t 2 ) = J d d xe 
where 

fo(x,h,t 2 ) 



1 ' " + ^ S -sin' 1 (f (x,t u t 2 )) (83) 



L T (ti, t 2 ) 



d d Q _ 



(2ir) d 



iQ-x 



(84) 



We will now show that F a (Q) falls off algebraically 
with large wavenumber. This means that the short- 
distance behavior is controlled, as in the OJK the- 
ory, by the auxiliary field m. We see that the cross- 
over functions H and Hq do not influence this short- 
distance behavior if we choose H = Hq. Notice also 
that it is advantageous to choose H(Q) to have a 
high-power of Q component so one can ignore the 
corrections to F(Q) = F a (Q). For this reason we 
keep the 710 term in Eq.(Q) and the crossover func- 
tion H(Q) does not influence the large Q behavior 
of F(Q) until terms of 0(Q~ 14 ). We do not expect 
significant variations in our numerical results to de- 
pend on whether we keep, for example, 710 or 73. 

Let us look in more detail at the large Q be- 
havior of F a (Q). This requires several steps. First 
we have, using Eq.(J52), that in coordinate space 



F a (x) = -sin l fo(x) 



(87) 



VI. LARGE Q AND SHORT-DISTANCES 

One of the main successes of the OJK theory 
is inclusion of the short distance nonanalytic behav- 
ior associated with Porod's law and the Tomita sum 
rule. We look here at how all of this fits into the 
present development. 

Let us consider Eq.(|5C|) for equal times ii = 
t 2 = t 3> to where 

F(Q) = F(Q,t,t) = ( ld J± e -V^ H ^ 
Jo Vi 



where fo(x) is given by Eq.(84) with ti = t 2 = t 

e -N(Q) 



fo(x) 



d d Q , 



(2ir) d 



Jd 



The short-scaled distance expansion for fo(x) can be 
written in the form 

f {x) = 1 - a x 2 (1 + b x 2 + c x 4 + d x 6 + • • •) . (89) 
If we define the integrals 

dQ Q n e~ N ^ (90) 



8 



then 



bo = - 



co 



14 
6 7 2 

2O/4 

7!J 4 



(91) 
(92) 
(93) 



F 2B = 47r(-l) 



n+l 



dx 2 



x=0 



F 6 = -96™/? 



-fs = 47r6!a 7 



(102) 



Using Eq. fl95|) in Eq.( |102| ) we easily find 

F 4 = 8na (103) 



(104) 



(105) 



9! J 4 



(94) 



Inserting Eq.(|89|) into Eq.(|37j) and expanding, it is 
only a matter of stamina to show: 

F a (x) = 1 - ax (1 + (3x 2 + 7a: 4 + vx & + ■ ■ •) (95) 

where 



a = — y2ao 

7T 



0: 



an 
12 



Co 



7 2 8 



60 
2 

a b 



3 2 
160 a ° 



(96) 



(97) 



(98) 



d o a 3 , 4 2 



25 a2 9 

-/3 a + - 7 a - 7/? 



(99) 



Thus all of the expansion coefficients for F a {x) are 
known in terms of the parameters A and fi in N(Q). 

The large-Q behavior of F a (Q) follows in 
three dimensions from the Fourier representation 

4tt Z 100 

Fa(Q) = 77 / dx x sin(Qx)F <T (x) (100) 
Q Jo 

and, after repeated integrations by parts, one ob- 
tains for large Q 



where the Porod coefficients are defined by: 



(101) 



Mo 



-4tt8!c 



(106) 



Following Tomita we can show that the lack 
of even terms in the expansion of F(x) and F a (x) 
leads to the set of sum rules 



Sn — 



d 3 Q 
(2tt) 3 



Q 2 



n+l 



F(Q) E 



F2m 
C\1m 



m=2 



= . (107) 



If we set F(Q) = F a {Q) in Eq.(|l07|), because 
of Eq. (|9^) , the sum rules are obeyed identically. 
However, if we insert F(Q) given by Eq.(^5|) into 
Eq.(107), there is no reason in general to expect the 
sum rules to be satisfied. Indeed this gives a set of 
conditions on F(Q) which can be used, along with 
the normalization 



So = 



d 3 Q 
(2tt) 3 



F(Q) -1 = 



(108) 



to determine the parameters A, fi and those deter- 
mining H(Q). 



VII. DETERMINATION OF THE SCALING 
FUNCTION 

The equal-time scaling function F(Q) given 
by Eq.(|50|) is a function of the 7„, 7°, A and fi. Our 
basic assumption is that our model is characterized 
by the normalization Eq.(lOS) and the set of sum 



rules given by Eq.(107). Thus we have an infinite 
set of parameters and an infinite set of conditions. 
Here we work out the truncated theory where we 
use the four conditions 5, = 0, for i = 0,1,2,3 to 
determine the parameters P — {72, 710, A, /x} with 
7„ = 7„- Thus there are no free parameters. 

One then has a rather complicated numerical 
minimization problem. First assume values of the 
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four parameters P^ and compute F^(Q) and si- 
multaneously the four sum rules Sn . Then choose 
another set P^ and determines the set Sn\ One 
then needs a measure, like 

^ W =E(^ l) ) 2 , (109) 

to minimize. Thus, if < J^ x \ the set of param- 
eters P^ is preferred over the set P^\ One then 
iterates this process until it converges to a selected 
fixed-point set of values of the parameters P* and 
scaling function F*{Q). 

In the numerical determination of the S n it is 
important to take into account that these integrals 
are slowly converging for large wavenumbers. Sup- 
pose we have integrated a sum rule out to a cutoff 
wavenumber Qm and obtain a contribution S$ M . If 
Qm is large enough, the system is dominated by the 
large-Q power-law behavior as in Porod's law and 
one can evaluate the remaining contribution to the 
integral in terms of an appropriate Porod coefficient: 

Sn=S^+^f- - (HO) 

The determination of the S® M is a challeng- 
ing numerical problem since one must perform mul- 
tiple nested integrations y\ , y 2 and internal Fourier 
transforms and still maintain sufficient numerical ac- 
curacy that the Porod coefficients F2 n can be ex- 
tracted and the sum rules constructed. Thus the 
choice Qm must be monitored carefully when deter- 
mining S^ M . 

As a result of extensive iterations we arrive 
at the fixed point values (3|] for the parameters: 
72 = 1.4994 . . ., 7i = 0.1232 . . ., A = 6.395 . . ., and 
fi = 0.2162 .... In Fig. 1 we plot the determined 
scaling function F(Q) versus the accurate numeri- 
cal results of Oono and Shinozaki ]22^ ]. Notice that 
the structure factor is normalized by the position 
and height of the first maximum. The agreement for 
small Q is very good and the width of the peak is in 
good agreement with the numerical results. In Fig. 2 
we plot Q 4 F(Q) for the theory and the same set of 
numerical results. Again overall agreement is good. 
The value of the Porod coefficient F4 for the the- 
ory and simulation are in excellent agreement. The 
feature of a second maximum is present in the the- 
ory but its position and width are not in very good 
agreement with the numerical results. In Figs. 3-6 



we plot the running contributions to the sum rules 
to give a feeling for the degree of convergence of 
the numerical procedure described above. To ob- 
tain the complete contribution to the sum rules we 
must add the last term in Eq. ( |llC| ) . The analytically 
determined Porod coefficients needed in completing 
the analysis are given by F4 — 21.87, Fq — 19.72, 
F s = -52.85 and F w = -25.41. 

The results here should be compared with 
previous direct calculations of the F(Q). In Ref.( 
p| ) the calculated F(Q) obeyed Porod's law, did 
not obey the lowest Tomita sum rule, went as Q 2 
for small Q and gave a width significantly too broad 
compared to numerical results. In Ref.( |l2|) the 
large and small Q limits were in agreement with ex- 
pectations but the Tomita sum rule was not satis- 
fied and the width of the scaling function was sig- 
nificantly too narrow. In Ref . ( f3lf| ) the authors ob- 
tained a scaling function which satisfied Porod's law, 
the Tomita sum rules and was fit to the numerically 
determined width. Unfortunately, this work not only 
did not give the Q 4 small Q behavior for F(Q), it 
did not satisfy the conservation law, F(0) > 0. 

VIII. AUTOCORRELATION FUNCTION 
EXPONENTS 

We turn next to a discussion of the two-time 
order parameter correlation function. In this case 
we focus on the behavior of the autocorrelation func- 
tion, which, when both t\ and ti are in the scaling 
regime, is given by 

*(*i,i 2 ) =< <Mr,iiMr,i 2 ) > 

2 f d d Q f^ V3 d Vl 3/2 . 

J ( 27T r J(to/T)V3 vi 

x — U q {t2,Tyl ,2 )H {Qy/y^Ho{Qy/m) 

J(t /T) 2 /3 J/2 

x£ d (y 1 ,y 2 )F a (Q£(y 1 ,y 2 ),Tyt /2 ,Ty 3 2 /2 ) . (Ill) 

In the regime t\ 3> t%, ^(ti^) is expected to take 
the form Q| 
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*(*l,t 2 ) 



L{ti) 



(112) 



and the exponent A is distinct from the growth law. 



It is not difficult to extract A from Eq.(lll). 

For ti ^S> t 2 the integral over y 2 is restricted 
to small values and the autocorrelation function can 
be put into the form 



*(ii,t a ) =V>o 



d d Q 
(2ir) d 



o J/i 



for the field at ii and the initial correlation. Then 
one has 



C(q,t 1 ,t ) 



di U q {h,l)N q {T) < a q {t)a- q {to) > (118) 



The analysis of this quantity for large t\ follows the 
earlier analysis with respect to rescaling of the t\ 
dependence. The main difference in the calculation 
is that there is no y 2 integration as in the t 2 3> to 
case and one has the direct factor as in Eq.(llt) 
which leads to the behavior 



x Ho{Q^yl) [ 
Jo 



2(i 2 /*l) 2/3 



— 72Q y-2 

V2 



xI d (y 1 ,0)F a (Qe(y 1 ,0),Ty 3 1 / \Ty 3 2 /2 ) . (113) 

They key observation is that for t\ 3> t 2 the aux- 
iliary field correlation function is small and we can 
linearized Eq.(|3|) relating the order parameter and 
the auxiliary field correlation function to obtain 



F (T (Q,ti,t 2 ) « -fo(Q,h,t 2 ) 

7T 



(114) 



where 



A(0 , (1 , (J)H ™>y ,2 £^ . (115) 



Lrp 



Id 



Then to leading order for t\ ^> t 2 

fo(Q,h,t 2 ) 



2L{t 2 )\ d/2 e -M Q3+ ^ Q ~ AQ ) 



Id 



(116) 



Inserting this result back into Eq.(113), we find that 
the overall time dependence is governed by the y 2 
integral given by 



2(L(t 2 )/L( tl )) 2 



dy 2 yT » (i(t 2 )/L(ti)) d/2+2 (117) 



and we can identify |33j A = d/2 + 2. This result cor- 
responds to the lower bound established by Yeung, 
Rao and Desai E^ . 

If t 2 = to then the analysis must be altered. 
The key point is that one must keep the first term 
in Eq.(p4[) and the dominant term for long times is 
the correlation between the second term in Eq.(|24 



*(ti,ta) « (Li) 



-d/2 



(119) 



In this case we again obtain that the exponent is 
given by the lower bound value Ao = d/2 found by 
Yeung, et al fj. 



IX. DISCUSSION 

In this paper we have developed a theory of 
the phase-ordering kinetics of the CH model which 
plays a role similar to the OJK theory for the NCOP 
case. The theory includes all of the desired features 
discussed in section 2 including the elusive Tomita 
sum rule. It is also in good numerical agreement 
with the best available numerical results. The the- 
ory, because of the nonlocal mapping onto the auxil- 
iary field, is much more difficult to treat analytically 
when compared to the OJK theory. Nevertheless it 
seems to work well. While we are able to extract the 
apparently mean field values for the nonequilibrium 
exponents A and Ao, it looks very difficult to include 
the nonlinear terms in the m field necessary to give 
the higher-order corrections for the exponents. 

There are several directions in which this 
work can be extended. One can, in principle, in- 
clude additional j n parameters in the model and 
satisfy higher order sum rules. This will be numeri- 
cally difficult and probably will not improve the scal- 
ing function significantly. While the detailed anal- 
ysis here has been for three dimensions, with only 
a modest amount of additional analytical work, the 
selection process can be applied to two dimensional 
systems. The most interesting new direction is to 
apply this same theory to off critical quenches. The 
expectation is that the scaled structure factor broad- 
ens significantly as one moves away from critical 
quenches. It will be interesting to see how well the 
theory describes the system as one moves toward the 
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coexistence curve and the regime where the Lifshitz- 
Slyosov- Wagner theory pm is applicable. 
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FIG. 1. Plot of scaled structure factor versus scaled wavenumber. Solid line represent the theory and the squares 
the numerical results from Ref.([22]) 
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FIG. 5. Plot of sum rule S$ M versus cutoff wavenumber 
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FIG. 6. Plot of sum rule M versus cutoff wavenumber 
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